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We use a stochastic approach to show how Taylor dispersion is affected by kinetic processes of 
adsorption and desorption onto surfaces. A general theory is developed, from which we derive 
explicitly the dispersion coefficients of canonical examples like Poiseuille flows in planar and cylin- 
drical geometries, both in constant and sinusoidal velocity fields. These results open the way for 
the measurement of adsorption and desorption rate constants using stationary flows and molecular 
sorting using the stochastic resonance of the adsorption and desorption processes with the oscillatory 
velocity field. 
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I. INTRODUCTION 

In presence of a uniaxial stationary laminar fluid flow, 
the diffusion of Brownian particles along the direction of 
the flow is enhanced by an amount proportional to the 
inverse of the molecular diffusion coefficient D b . This ef- 
fect, known as Taylor dispersion [TJ[2], originates from the 
combination of the inhomogeneous velocity field experi- 
enced by the particles and the diffusive transverse motion 
which leads to a random sampling of these fast and slow 
streamlines. Taylor dispersion has implications in many 
fields, especially those involving chemical reactions and 
determination of reaction rates such as microfluidics and 
chromatography, and has thus been the subject of a num- 
ber of works both at the theoretical and experimental 
level [3-8 . In numerous practical situations, Taylor dis- 
persion in the bulk flow is coupled to the adsorption and 
desorption processes taking place at the walls confining 
the fluid. So far, the theoretical analysis of the resulting 
process has mainly been done explicitly in two limiting 
situations. 

In the first class of models, the transverse motion is 
not explicitly considered, which physically corresponds 
to the infinitely well stirred limit of high diffusion coef- 
ficient Dij. A representative example is the famous two- 
state model of chromatography introduced by Giddings 
and Eyring in 1955, in which a particle can be either in 
the mobile phase (in the flow) or in the immobile phase 
(adsorbed on the confining walls), the rates of change be- 
tween phases being constant [9] . An important extension 
concerns the case when the velocity of the mobile phase 
oscillates with time according to vcos(ut) [TUJ |TT] . In 
particular, stochastic resonance has been shown to oc- 
cur if the rates of change between phases are both equal 
to gj/2, leading to a maximum of the dispersion coeffi- 
cient [I2j [13]. This effect has recently proved to have 
applications in molecular sorting p~4] [15] . 
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The second class of models has investigated explicitly 
the transverse motion, but for specific kinetics of adsorp- 
tion and desorption: In [16] , the dispersion coefficient is 
calculated when the exchanges with the surface are in- 
finitely fast (local chemical equilibrium), while Biswas 
and Sen have considered the situation of irreversible ad- 
sorption on the surface [17]. Besides, these studies focus 
on stationary velocity fields and the important case of 
oscillating velocity fields mentioned above is not consid- 
ered. 

In this article, we develop a theoretical analysis of Tay- 
lor dispersion in presence of general adsorption and des- 
orption processes. Relying on a stochastic approach (i) 
we derive explicit expressions of the dispersion coefficient 
for the canonical examples of Poiseuille flows in planar 
and cylindrical geometries, both for stationary and oscil- 
lating velocity fields, thus opening the way to the deter- 
mination of heterogeneous rate constants from the mean 
velocity and dispersion coefficient; (ii) we recover the fact 
that, in the case of a stationary velocity field, the sources 
of dispersion associated to bulk transport and adsorption 
and desorption processes combine additively p~8l[T9]; (iii) 
in the case of an oscillatory velocity field, we show that 
the dispersion coefficient can be optimized and discuss 
possible implications in the context of molecular sorting. 



II. THE MODEL 

We consider a Brownian particle in a flow of velocity 
field v in direction x. The position of the particle in the 
transverse direction is denoted by y and the full position 
by r = (x,y). The longitudinal dynamics of the particle 
is assumed to be given by the Langevin equation: 

x(t) = v(y(t),t) + l b (y(t))r] b (t) + l s (y(t)) Vs (t), (1) 

where lb(y(t)) stands for the indicator function of the 
bulk b (equal to 1 if the particle's position is in the bulk 
and otherwise) which accounts for bulk diffusion (with 
diffusion coefficient D b ) and l s (y(t)) for the indicator 
function of the surface s associated to surface diffusion 
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(with diffusion coefficient D s ). The independent Gaus- 
sian white noises r] b and r] s are defined by their correlation 
functions: 



(Vb(t)) = (Vb(t)) = 
{Vb(t)Vb(t')) = 2D b 6(t-t') 
(f ls (t)r ls (t')) = 2D s 5(t-t'). 



(2) 



The transverse diffusion equation is driven by the evolu- 
tion equation: 

' d t P(y,t\y',0) = D b V 2 P(y,t\y',0),Vy&b 
d t T(y,t\y',0) = -k d T(y,t\y',0) + k a P(y,t\y',0) 
= D b d n P(y,t\y',0),\/yes 

(3) 

where k a (resp. kd) is the adsorption (desorption) rate 
in length- time -1 (resp. time -1 ) [20] , P (resp. V) is the 
propagator corresponding to a final state in the bulk 
(resp. on the surface) and d n stands for the normal 
derivative. Initially, the particle is assumed to start from 
x = and the process y(t) to be stationary, characterized 
by the stationary distribution P s tat(y) (uniform within 
each phase and depending only on the ratio k a /kdL) and 
the transition probability P(y, £|y', t') = P(y, t — t f \y\ 0). 

The first two moments of x(t) are then respectively 
found from Eqs. 0, ([2| to be given by 



(x(t)) 



f dt' fdyP stat (y)v(y,t') (4) 

JO Jb 



where the integration domain of the spatial integral is 
the transverse cross section in the bulk b and 

{x 2 (t)) = fdt' f dt"(v(y(t'),t')v(y(t"),t")) 
Jo Jo 

+ 2D b (T b (t))+2D 8 (T a (t)), (5) 

where T b (t) (resp. T s (t)) is the cumulative time spent in 
the bulk (on the surface) up to time t. Finally, one has 

(x 2 (t)) - (x(t)) 2 - 2D b (T b (t)) - 2D s (T s (t)) = 
2 [ dr [ dt f [ d Yl [ dy 2 v{y 1 ,t')v{y 2 ,t' - t) x 

Jo Jt Jb Jb 

x P s tat(Y2) [P(yi,r|y 2 ,0) - Pstatfrl)] = f(t). (6) 

We now specialize this general formula to the two 
experimentally relevant cases of either a stationary or 
sinusoidal velocity field. 



III. CASE OF A STATIONARY VELOCITY 
FIELD 

In the case of a stationary velocity field v(y, t) = v(y), 
the large time limit of the variance of the displacement 
is easily shown from Eq. $6§ to be given by 



where the dispersion coefficient is K = P stat (b)D b + 
P s tat(s)D s + K v , where P 3ta t{b) (resp. P 3ta t(s)) is the 
stationary probability to be in the bulk (resp. to be ad- 
sorbed on the surface) and the velocity-dependent part 
reads: 

K v = dyi / dy 2 v(y 1 )v(y 2 )P8tat(y2)h(y 1 \y2), (8) 

Jb Jb 

with 

POO 

ft(yi|Y2)= / [P(yi,t\y 2 ,Q)-P stat ( yi )]dt. (9) 
Jo 

Note that /^(yi|y2) is the pseudo-Green func- 
tion [21 of the transverse problem, which satisfies 
-A>V 2 /i(yi,y 2 ) = S(y 1 - y 2 ) - P s tat(yi)- As soon as 
this pseudo-Green function can be determined, Eq. (|8| 
provides a general expression of the Taylor dispersion 
coefficient in presence of adsorption and desorption 
processes. 

Importantly, this expression can be made fully ex- 
plicit in the canonical examples of planar and cylindrical 
Poiseuille flows, corresponding respectively to velocity 
fields v(y) = 6v^ (l — ^) (the transverse cross section 
being a segment of length L and y G [0, L]) and v{r, 0) = 

v{r) = 2v ^1 — -^2^ (the transverse cross section being a 

disk of radius R and y = (r, 0) e [0,R] x [0,27r]), where 
v stands for the velocity averaged over a cross section. 
The explicit determination of the pseudo-Green func- 
tion h(yi\y 2 ) is conveniently performed by first Laplace 
transforming Eq. (|3| , calculating the Laplace transform of 
the propagator and then going to the small Laplace vari- 
able limit. Note in particular that the boundary condi- 
tions Eq.(|3| associated to adsorption/desorption become 
simple radiative boundary conditions in the Laplace do- 
main (see for example ref. [22 J. Lengthy but straightfor- 
ward calculations finally lead to symmetrical functions of 
their arguments, which read 



h pl ™(yi\y 2 ) = 



^? + feyi + ^ 2 2 -y 2 (i + fe) 



1 r3 , 2D b k a 

3 L +-fcT" 



k a L z 



d > ( L + 2 -ty 



(10) 



if Vi < y-2 and 



1 -In-* 



— / de 2 h^-(r 1 ,e 1 \r 2 ,0 2 ) 

27r J 2nD b r 2 

(2R + 4^) (r? + rf ) - 3R 3 + 16^ 



8ttRDj 



(x 2 (t))-(x(t)} 



2 ~ 2Kt, 

t— >-00 



(7) 



(ii) 
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if ri < r2 . The velocity- dependent part of the dispersion 
coefficients is then found to have the same form in both 
geometries: 



K 



Pois. 



a- 



ka 



kd 



(12) 



where the length / and constants (a, /?, 7) are to be sub- 
stituted by L and (^0,102, 18) in the planar case and R 
and (^,44,12) in the cylindrical one. 

A few comments are in order, (i) The specific case of 
infinitely fast exchange with the surface (local chemical 
equilibrium), considered in the cylindrical case in [16] as 
one of the generalizations of Taylor dispersion, is recov- 
ered in the joint limit k a —> 00, kd — » 00 with k a /kdL 
fixed, and is given by the first term of the right hand 
side (r.h.s.) of the general expression (12). (ii) In the 
infinitely well stirred limit D b — > 00, K^ ols - is reduced to 
the second term of the r.h.s. of Eq. (12). This second 



source of dispersion, associated to adsorption and des- 
orption kinetics only, corresponds to the 0-dimensional 
dispersion coefficient used in the usual two-state model 
of chromatography [12] , with a rate 2k a /l (I = L,R) for 
the transition from mobile to immobile states, (iii) Note 
that these two contributions to dispersion turn out to 
combine themselves additively. (iv) The knowledge of the 
mean displacement (x(t)) ~ P stat (b)vt (see Eq. Q) and 
the dispersion coefficient Eq. ( [l2] ) enables the experimen- 
tal determination of k a and kd from the measurement of 
the average velocity and dispersion coefficient for a given 
species interacting with a given surface, (v) It also allows 
to design a chromatographic column for the separation of 
a mixture with known adsorption and desorption rates. 
However, in the important case where the components 
of the mixture have different kinetic rates k a and kd but 



similar partitioning coefficients k a /kdL, such a separation 
is impossible with a stationary flow, since P s tat(b), hence 
(#(£)), of all species are similar. As we now proceed to 
show, efficient sorting can nevertheless be achieved in this 
case by resorting to an oscillatory flow and exploiting the 
stochastic resonance between the exchange kinetics and 
the flow, thereby extending the idea put forward by Alcor 
et al. in another context [14[ [15] . 



IV. CASE OF A SINUSOIDAL VELOCITY 
FIELD 

We now consider the case where the velocity field is a 
sinusoidal function of time v(y, t) = v(y) cos(ujt). In this 
case, the average position (x(t)V tends to at long time, 
while the function f(t) in Eq. (|6| becomes: 

/(*)= f dr [ dt' I dy x [ dy 2 v(y 1 )v(y 2 )P s tat(y2) x 

JO Jt Jb Jb 

x [^(yi,^|y2,0) -P s tat(yi)} [coscj(2t' -r) + cos ur] . 

(13) 

In the large time limit, the term in coso;(2t / — r) be- 
comes negligible and the velocity-dependent part of the 
dispersion coefficient is given by 



Kv = \ J b dyi J b dy2V ^^ v ^ 



lb Jb 

xPstat(y2)Re(P(yi, -Hy 2 )), (14) 

where Re(P(y±, — iu\y2)) stands for the real part of the 
Laplace transform of the propagator, with the Laplace 
variable s = — iuo. We focus here on the canonical case 
of the planar Poiseuille flow described above. The final 
result for the velocity-dependent part of the dispersion 
coefficient reads: 



K 



Pois 



_ L 2 v 2 3Z xt cosn ^ + X c cos X + \f smn X + x s sm X 



D b 2X 6 (2Y + Z) pt cosh X + P c cos X + pf sinh X 



where we have introduced the polynomials 



- p s sin X 



(15) 



xf = 2X A Y 2 - 6X 2 Y ± (X 2 Z 2 + X 2 - 12Y) 
X f = 2X 3 YZ + 12XY ± (2X 3 Y - 12XY 2 - 3XZ 2 
pf = 2X 2 Y 2 ± (Z 2 + 1) 
pf = 2XYZ ± 2XY 



\2XYZ - 3X) 



(16) 



and the reduced variables X = L^/u/2Di ) ^ Y = k a /ujL 
and Z = kd/oo. This expression constitutes one of the 
main results of the present work. It extends in particular 
the results known in absence of adsorption and desorp- 
tion at the wall [12] (recovered in the limits Y — >• or 
Z — >• 00) and in the case of an infinitely well stirred limit 



in the transverse direction (corresponding to X — >• 0) 
HQ] [H[T4] [TS , with an adsorption rate 2k a /L. 

Besides being an important theoretical result in itself, 
the expression (15) also allows one to discuss the possi- 



bility of sorting components of a mixture in an oscillatory 
flow, in presence of adsorption and desorption processes. 
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Alcor et al [HI [15] have demonstrated, both theoretically 
and experimentally, that an oscillatory driving (by an 
electric field) can be exploited to separate components of 
a mixture that switch between two bulk states with differ- 
ent rates. Using the O-dimensional two-state model gen- 
erally considered in chromatography, they showed that 
the dispersion coefficient displays in this case a maxi- 
mum when the two rates are close to cj/2, i.e. when the 
average time spent in each state are equal and compara- 
ble to a half-period of the driving, resulting in an effective 
rectification of the flow experienced by the particles. We 
now discuss the possibility of extending this idea to the 
case of Taylor dispersion with adsorption and desorption. 



mum K° pt . The optimal rates satisfy: 



Log, [Z] 




Log, [Z] 



LogiotZ] 



FIG. 1. (Color online) Velocity-dependent part K v of the 
dispersion coefficient, normalized by the value at the optimum 
Ky pt as a function of the reduced rates Y = k a /ooL and Z = 
kd/uj, for three values of X = L^/uj /2Db'. 0, 1 and 3.5. 

Importantly, we find that Taylor dispersion with ad- 
sorption and desorption at the walls of the confining 
surfaces also exhibits a stochastic resonance. However, 
the picture is more complex due to the coupling between 
motion in the direction transverse to the flow and the 
exchange processes. Indeed, K^ om \uo) displays a global 
maximum for a set of finite reaction rates (k a ,kd), pro- 
vided that X < X* w 3.5, as can be seen in Figure [T] 
which reports K v normalized by the value at the opti- 



k° pt =ujLY opt (X) and k 



opt 



jZ opt (X) (17) 



where y opt and Z opt , reported in Figure [2^i, are well 
approximated by their small X expansions for an appre- 
ciable range (see Figure 



^opt 



and 



11 

336' 
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50400' 



X 4 



331 



1241856 



-X" 



G(X 8 ) 
(18) 



jjopt _ ± 



71 

840' 



17 
1800' 



-X 4 



25427 
25872000 



X 6 + 0(X 8 



(19) 

and correspond to the following expansion of the optimal 
dispersion coefficient: 



K°*\X) 



L 2 v 2 
D h 
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79X 2 



32X 2 3360 564480 



0(X 4 ) 



(20) 




FIG. 2. (Color online) (a) Optimal reduced rates, y opt = 
k° a pt /ojL and Z opt = k° d pt /oj as a function of X = L^u/2D h 
maximizing the velocity-dependent part of the dispersion co- 
efficient K v . (b) Corresponding values of the optimum if° pt , 
in units of L 2 v 2 /Db- Numerical results (solid lines) are com 



pared to the analytical expressions Eqs. (18)-(20) for small X 
(dashed lines). In (b), K v in the absence of adsorption and 
desorption is also indicated (dotted line). 

While for small values of X, the resonance is obtained 
for kd ~ 2k a /L ~ cu/2 and has similarities with the case 
of Alcor et al. [TU [15] , the optimal rates can in fact dif- 
fer by orders of magnitude from uj as X — >• X*. More- 
over, while the resonance process is rather selective in the 
small X regime, as K v decreases rapidly around the opti- 
mum, Figure [l] indicates that the selectivity deteriorates 
when X increases. Finally, the optimal dispersion coef- 
ficient, enhanced by orders of magnitude by adsorption 
and desorption at the walls, also decreases with increas- 
ing flow frequency (see Figure ^p). All these observations 
reflect the transition from the well-stirred regime X — » 



5 



where "chemical" and Taylor dispersions combine addi- 
tively (the first two terms in Eq. ( 20 ) are proportional to 
v 2 /k^ pt and v 2 L 2 /D^,, respectively), as in the stationary 
case, to one where exchange at the walls can be limited 
by the diffusive influx of mobile species. In the latter 
case, only a fraction of mobile species can benefit from 
the rectification mechanism allowed by the adsorption 
and desorption processes. Its effect on dispersion is then 
maximal when the particles spend the same time at the 
surface and in this boundary layer accessible by diffusion 
within a flow period (now regardless of reaction rates 
provided that the exchange is fast compared to the flow 
frequency) . 



V. DISCUSSION 

This optimization of the dispersion coefficient can be 
used in the context of molecular sorting, where the prob- 
lem is to isolate a species of given (fc a , kd) from a mixture. 
In fact, the ratio Z°^/Y°^ = k° d pt L/k° a pt = g(X) is an 
increasing function of X, such that 2 = 7_ < g(X) < 
7 + ~ 5 for X < X*, so that the parameter L can be cho- 
sen close to 2k a /kd. This ensures that X = g~ 1 (Lkd/k a ) 
is arbitrarily small and determines the corresponding 
value of uj = 2DbX 2 /L 2 to tune. Finally, this shows that, 
in principle, one can always determine parameters L and 
uj such that K^ ols - is maximal for (k a ,kd) while main- 
taining a small value of X to ensure efficient sorting, and 
thus define an optimal setup to be used. Furthermore, 
our approach quantifies the theoretical efficiency of such 
sorting. 

In practice, experimental constraints limit the acces- 
sible range of system sizes (L m i n , L max ) and flow fre- 
quencies (cj m i n , cj max ). As an example, for a typical 
microfluidics setup, these ranges are of the order of 
(10 /im, 1 mm) and (0 s _1 ,10 s _1 ), respectively. Sev- 
eral cases have then to be considered: (i) For species 
such that L m i n kd/k a < 7_ < L mdL ^kd/k a ^ the optimal 
setup for efficient sorting defined above is indeed realiz- 
able and the constraint on uj is irrelevant for all values of 
(ii) For species such that 7_ < L max kd/k a < 7 + , the 
maximum sorting efficiency can be reached for L = L m [ n . 
This implies a corresponding value of uj = ID^X 2 / / ^ min , 
which constraints the range of applicable diffusion coeffi- 
cients Db to Db < uj rnax L 2 nin /2X 2 . (hi) For species such 
that L min k d /k a > 7+ or L max k d /k a < 7_, the method is 
in practice not applicable. 

Microfluidic technology offers a particularly versatile 
set of tools to taylor the geometry, flows and surface prop- 
erties to design experimental setups for the measurement 
of adsorption and desorption rates (stationary case) or 
separative applications (stationary and oscillatory cases) 



according to the predictions of the present work. In prac- 
tice, for a typical microfluidic channel one has a trans- 
verse length L ~ 10 -4 — 10 -3 m and velocities up to 
v ~ 10 -4 — 10 -3 m.s -1 , while colloids, macromolecules 
and molecular solutes have diffusion coefficients in the 
range ~ 10 -12 — 10 -9 m 2 .s _1 . Oscillatory flows can be 
considered of the form v(y) cos ujt if momentum diffusion 
in the direction transverse to the flow is fast compared 
to the period of the flow, i.e. uj <C L 2 jv with v the kine- 
matic viscosity of the fluid. Such a condition is always 
satisfied for water in a microfluidic channel, given that 
^H20 °° 10 -6 m 2 .s _1 and that only frequencies smaller 
than uj/2tt < 1 s _1 can be achieved. This last point also 
indicates that the present approach will allow to measure 
sorption/reaction rates slower than Is -1 . As an exam- 
ple, the dissociation rate of DNA double-strands, which 
depends on the number of base pairs (bp), can be in the 
range 10~ 5 - 10" 3 s" 1 for a few tens of bp [23H25]. If 
one considers a surface grafted with single-strand DNA, 
one could selectively separate from a solution a strand 
containing the complementary sequence by adjusting the 
flow period and the grafting density in order to tune the 
adsorption rate k a (in the low surface coverage limit, it 
will be proportional to the latter). 



VI. CONCLUSION 

In conclusion, the present study introduces general 
analytical results which extend previous works on Tay- 
lor diffusion without adsorption and desorption or on 
a 0-dimensional two-state model valid only in the per- 
fectly stirred limit. This approach is not limited to the 
Poiseuille flows considered here as an illustration. In par- 
ticular, it can be straightforwardly extended to the case 
of electro-osmotic flows, in which an additional length 
scale, the Debye screening length, can be tuned by chang- 
ing the ionic strength of the solution. The now well- 
established microfluidic and rising nanofluidic technolo- 
gies offer a particularly versatile set of tools to taylor the 
geometry, flows and surface properties to design exper- 
imental setups for the measurement of adsorption and 
desorption rates (stationary case) or separative applica- 
tions (stationary and oscillatory cases) according to the 
predictions of the present work. 
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